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ABSTRACT 


AN INVESTIGATION OF THE SIDE THRUST AVAILABLE FROM 
PROPELLER SHROUDING SURFACES 


by 
SAMUEL JAMES GORDON 


Submitted to the Department of Mechanical Engineering and to the Department 
of Naval Architecture and Marine Engineering in partial fulfillment of the 
requirements for the degree of Master of Science in Mechanical Engineering 
and the Professional Degree of Naval Engineer. , 


À theory and analytical method are developed for the determination of the 
circumferential circulation distribution, over an annular airfoil of arbitrary 
included angle. The resultant circulation distribution is then used to determine 
the forces acting on the airfoil surface. 


The theory is based on the "L-Method" developed by J. Weissinger, which 
consists of replacing the lifting surface by a single concentrated ring vortex and 
an associated system of trailing vortices. The method presented here is com- 
pared to other theories in the limiting case of the complete symmetrical annulus. 
It is shown to give the same result as Weissinger found using the "L-Method. " 
This, in turn, is shown to be a very close approximation to the more accurate 
theories and to experimental results. 


A sample calculation is made, based on the method presented, of the forces 
acting on a semicircular propeller shroud, to determine its usefulness as a 
steering device in lieu of the conventional ship's rudder. The results are com- 
pared with those for a fully movable spade type rudder with area equal to that 
of the shroud and with a high lift to drag ratio. It is found that the predicted side 
thrust from the shroud is more than twice the maximum predicted for the rudder. 
Additionally the shroud is found to give a thrust augmentation equal to approxi- 
mately 10% of the propeller thrust while the rudder absorbs a minimum of 1% 
of the propeller thrust and a maximum of 30%. 


Thesis Supervisors: Justin E. Kerwin 
Associate Professor of 
Naval Architecture 


Hal L. Moses | 
Assistant Professor of 
Mechanical Engineering 
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Introduction: 
The literature contains a good deal of material on the theory of annular 


(2) 


airfoils and ducted propellers. Authors such as Morgan and Weissinger(2 
have treated the problem of a symmetrical ring airfoil at an angle attack, 

alone, as well as in the presence of other objects such as propellers and hubs. 
Weissinger has also treated the case of a complete annular airfoil with deflected 
control surfaces. There seems to be no evidence of any previous work on in- 
complete rings. 

Since the incomplete shroud seemed, from an intuitive standpoint, to 
provide a very effective means of steering a ship, I constructed a small self- 
propelled model and fitted it with a semi-circular propeller shroud. The 
model is illustrated in Figure 1. By rotating the shroud about the propeller 
axis the side thrust developed by the shroud could be directed so as to cause 
the model to turn to the right or left. The steering device worked very well on 
the model and provided the stimulus for this investigation. 

The theory developed here, I feel, provides a means for making a reasonably 
accurate prediction of the performance of either complete or incomplete shrouds 


with a circumferential variation in circulation distribution. 





Conclusions: 

The results of this sample calculation indicate that partial shrouds or 
complete shrouds with controllable circulation and lift distributions can be used 
to good advantage as steering devices for fluid born vehicles such as ships and 
ground effect machines. 

The following improvements in performance over that obtained with con- 
ventional rudders are cred 

l. Substantial elimination of rudder drag at all turning rates. 

2. Thrust augmentation from the shroud. 

3. Improved steering capability with little or no way on. 

4. Significant decrease in turning diameter with an equivalent size 

turning device. 

The following possible shortcomings in the application of shrouds as 
steering devices require consideration and investigation: 

I. Decrease in propeller loading necessitated by shroud induced 

cavitation. 

2. Propeller induced shroud vibration. 

3. Possibility of loss of steering capability with engines stopped. 

It is hoped that further theoretical work in this area along with experi- 
mental testing will be carried out and that the apparent benefits to be gained 


from the use of propeller shrouds for steering will be employed in practice. 








SIDE VIEW OF MODEL 





VIEW SHOWING SHROUD IN POSITION FOR TURNING 
TO THE RIGHT 


FIGURE I 





Procedure: The method utilized for determining the forces on the shroud is that 
developed by J. Weissinger and known as the "L-Method" or "Weissinger's 
Method, |“ The following basic assumptions are made: 

(a) The fluid is incompressible and baratropic. 

(b) Body forces, such as gravity, are neglected. 

(c) The flow is steady. 

(d) Viscous effects are confined to a negligibly thin boundary layer 

which does not separate on the shroud. 

The method consists of the replacement of the lifting surface by a single, con- 
centrated bound vortex at the 1/4 chord point and an associated system of 
trailing vortices as shown in Fig. 2. The circulation distribution, along the 
bound vortex, is determined by requiring that the velocity normal to the mean 
line atthe 3/4 chord point must be zero. Since the circulation distribution varies 
along the span of the shroud, a system of trailing vortices will be shed off which 
will form a cylindrical sheet of radius equal to that of the concentrated vortex | 
ring. The boundary condition is linearized by requiring that the quotient of the 
radial velocity and the free stream velocity, at the 3/4 chord point, along this 


cylindrical sheet, be equal to the slope of the section at the 3/4 point. That is: 


dr. az dr, 
—— "UN = 712,6, 13/4) -a (9) | 


"BO Q5 Aga) 


Once the circulation distribution is known, the radial force and axial forces 
per unit length at any section can be found by the "Kutta-Joukowski Law" which, 


in this case, states: 
P2: F,(8) = "P Vyp(l, 0, 0) L(g) 


P3: F,(6) = elar, (1, O, 0) + er (1, O, 1/4) Tie) 
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BOUND VORTEX 
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TRAILING VORTEX 


ze 


(Thickness not shown) 
FIGURE 2 


Shroud length. 

Limiting half angle of shroud. 

Angular coordinate along concentrated vortex in 1/4 chord plane. 
Radius of concentrated vortex in 1/4 chord plane. 

Pivotal point at radius equal to rj 4 in 3/4 chord plane. 

Angular coordinate of pivotal point in 3/4 chord plane. 
Displacement of mean line relative to nose tail line, nondimensionalized 
on L. 

Slope of mean line relative to nose tail line. 

x/r]/4 = Nondimensionalized axial coordinate 

r/r) /4 = Nondimensionalized radial coordinate 

Convergence angles of nose-tail cone 


Angle of attack in plane of symmetry. 








For a shroud which is symmetrical about 8 = O, the resultant radial force 
R lies in the plane of symmetry and is equal to the total lift acting on the shroud. 
R can be found by integrating the components of F,(9), parallel to the plane of 
symmetry, over the span as follows: 


a 
Pg: T(0) = Vy jg J, Fre) cos 0 do 


Similarly, the total induced drag D, acting on the shroud, is found as 


follows: 
De (3) 
PS: Dis Vy jg J Fx(0) «8 


The problem reduces to the determination of the circulation distribution 
for a given propeller geometry and loading and a chosen shroud geometry. 


The circulation distribution is represented as a Fourier cosine series. 


oo 


A-23: I (g) m ) Án cos (n0) 
n=0 


It is to be noted here that 8 is the actual polar angle measured from the plane 
of symmetry. For an incomplete shroud the condition that the circulation 
vanish at the tips is satisfied by requiring that the sum of all terms in the series 
vanish rather than that each term vanish individually. “The resultant radial 
velocity induced by the shroud, at a point on the trailing vortex sheet, is deter- 


mined as a series equation in terms of a set of coefficients and integrals. 


6 
A-26: de. = 7) A 


n=0 m=l 





10: 


where Cui are coefficients involving the parameters Ap and L m are geometric 
integrals. 

Since four terms in the series representation of l'(g) (Eq. A-23) will 
provide a representation of more than sufficient accuracy for our purposes, we 
will restrict our consideration ton 2 0, ..., 3. We must, therefore, evaluate a 
maximum of 24 coefficients, C,,,, and 24 integrals, L,,, for each angular co- 
ordinate (¢) at which we desire to determine dre" These can then be combined, 
using Eq. A-26, to form an equation in terms of the coefficients, An, of Eq. 
A-2839. If weletn=0,..., ee SU DES 3, we must evaluate dr. at (nmax 1) 
values of $ by using Eq. A-26. 

Knowing the propeller velocity field for a given propeller geometry and 
loading, and knowing the shroud geometry we can evaluate Ar; at the same 
(Mmax 7 1) values of $ by Eq. P-1. In the case of an incomplete shroud we have 
an additional boundary condition to satisfy and that is that I. £8) ^ 0. This gives 
us an additional equation resulting in a set of nmax simultaneous equations in the 
nmax Unknowns Ag» +++» An max: For a complete shroud we must evaluate 
dr. at an additional point to give us nmax equations. 

The solution to this set of equations provides the desired representation of 
the circulation distribution and permits the determination of the side thrust 
available from the propeller shroud as well as the drag induced by it. 

Weissinger has shown,in Ref. 1, and it is shown here in Appendix C, that 
this method yields the exact value of total lift for a complete ring airfoil in the 
limit of very high or very low aspect ratios (L/2rj/4) and yields a remarkably 
good approximation for arbitrary aspect ratios falling between these limits. 

It is to be noted that this method gives the correct solution only for airfoils 
of zero thickness as it takes account only of the circulation and not of the thick- 


ness. The effects of thickness must be accounted for separately; however, their 





e 


effect on the total force acting on the surface is small in comparison to that of 


circulation as shown in Fig. 22, Ref. 2. (Reproduced here as Fig. A-5.) 





EE 


Sample Calculations: 


Consider a propeller with the following characteristics: 


Dp = 12.0 SHP = 6000 
(1-t) = 0.80 (1-W) = 0.70 


Vs = 20 KTS V, = Vo 5 14 KTS; V, = 23.65 F. P. S. 
My = 0.625 epr = 1.00 
N = 200 RPM 
" u — TVa 
THP - SHP —nge,, - 3750 HP = ——- 
Cp = Thrust Coeff. = B 
DET 2 
TRp 2 ev, 
Cr = 2 THP) x550 ; p21.99 slugs/ft? 
T Rp e(V..) 
2(3750) x 550 T 7.8 x 8.5 x 109 


E 
m(6. 0) (1. 99) (23.65)> m(3.6) (1. 99) (2. 365)? x 10% 


Cr = 1.385 


Assume the circulation distribution over the propeller is constant; radial and 


axial velocity fields for this case are tabulated in Tables 3 and 4, Ref. (8). 





E 


Shroud Geometry: 


Consider a shroud of constant chord, camber and convergence angle, and 
assume, for the purposes of this calculation, that the velocities induced by the 
propeller, at the shroud vortex sheet, are equivalent to those induced on a 


cylinder of radius Rp. Take: 


GE) = 0.85 L = 9.6" 


Choosing NACA a = 0. 6 meanline as tabulated in Ref. (7) 


= = + 0. 1827 


A = 2.58% ; c| x = 


L. 
rig 0. 0737 


|x= 


mie! 


Referring to Table 3, Ref. (8) at ü =1.0; == -0.4 (shroud leading edge) 


X 
'R 
— 


a = +0.08 ; d - 0.08 (1.385) v. 


a +Y =1.11y 


Referring to Table 4, Ref. (8): 


— 


E == 0.086; im - -0. 086 (1. 385) v. 


Vo CT 


Ne 








34: 





E 

VreTAL NOSE-TAIL TRAILING VORTEX 

-frp \ LINE SHEET 

dut vs = Es 

| I. MEAN LINE 
r « 
la ý 
- == =æ - ee xX 
FIGURE 3 


For ideal angle of attack: 


ap = -B+oj = an c zt CHE 


a Bta = 6. 12° 
-B = 6.12% - 2.58% = 3. 54° = 0. 620 Radians 


The shroud will have the following characteristics: 


rj/4 7 6.0’ + (0.737) (9. 6) + 0. 62 (9.6/4) = 6. 857% 


L = 9.60’ 
Q9 = 90° 
Mean line: N.A.C.A. a=0.6 
-8 = 3,540 
Br 4. 80 


^3/4 7 rij4 ^ 6.857 — 9 0 





"Sr 
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Boundary Conditions for Determining Thrust: 


(1) ato = + 9 T'(g) = 0 


de. +4 , 
2) at 374? ~R =C ^ c (-B)= 0.1827 - 0.062 
Up + Vo 
Des (xp + Vo) (0.1207) - dip 


From Table 3, Ref. (8): 


q 
eP = 0,155 
Veo Grp 


dyp 7 0.155 (1. 385) v,, = 0.215 v,, 


From Table 4, Ref. (8): 


ED. -0. 136 
Vo Cr 
Ip ” -0.136 (1.385) v_ 


dug = (0.147+ 0.188) v,, = 0. 335 v,, 


q, = 0.335 (23.65) = 7.93 F. P.S. 


Take: 
$7 = 0 


$) = 11/4 


$3 = 11/3 





-16- 


The set of simultaneous equations determining the four coefficients in the ' 


cosine series representation of the circulation distribution are: 


0.025 Ag + 0.0741 A, + 0.0963 Ay + 0.1739 Az = 7. 93 
0.0235 Ay + 0.0422 A + 0.0198 A, - 0.0450 Ag = 7. 93 


0. 0214 Ag + 0. 0248 A, + 0.0165 Ay - 0.153 Ag = 7, 93 


1 
l0Ag +0A¡  —10A3 +0Ag  =0 


Solving simultaneously: 


Ag = Ay = 338.0 
Ay = -249,2 
Ag = -84.6 
l'(9) = 338 - 249.2 cos (9) + 338 cos (29) - 84.6 cos (39) 


By the Kutta -Joukowski Law 





By Ref. (8): 
Nem = 0.10 ; Ip = 0.1385 Va 
X=O 
LoT) = -188, 500 Lbs. 


LoT) is positive when directed radially outward, perpendicular to Vi Sea 


NP 
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Determination of the Induced Drag: 
D¡(aT) = - T(0) tan a, 


D¡(aT) = + 188,500 (0.045) = 8,480 Lbs. 


The induced drag is in a direction parallel to Viota] . 





FIGURE 4 


-T(0) = "Lot COS Am + D, (aT) sin aT 
-T(0) = 188,500 (0. 9943) + 8480 (0. 1067) 
-[(0) = 187,500 + 940 = 188, 440 

E. “Ty © Leer) = 188,500 Lbs. 

-D; (0) = “Leary Sin (aT) = D. (cT) cos aT 


-D;(0) = 20,100 - 8, 400 = 11, 700 Lbs. 








Se 


Form Drag: 
For an N. A. C. A, 66-210 section 


By Ref. (7): 


Cp (Y a=0 = 0.0094 
where 


D 


Cp| = 
aso l (V +q y nr IE 
3 co xp 1/4 


2. D¿ = 0. 0094 (1. 0) [(1. 215) (23. 65) 1 (6. 85 7) (9. 6) 


D; = 9.4 (0. 826) n (6. 857) (9. 6) 
De = 1600 Lbs. 


Taking Df as acting in the +x direction 


Final Predicted Thrust from this Shroud: 
Ahead thrust = 9,000 Lbs. 


Side thrust = 188,000 Lbs. 





E 


Comparison with a Conventional Rudder: 
Assuming that a/6 = 5/7 where: 


a= angle of attack 


ô = rudder angle. 


For a rudder of Aspect Ratio = 3,for the best L/D ratio and of area equal to the 


area of the shroud: 
A = area = (6.86 + 0.48 ) 7 (9. 6) = 221 Ft. 2 
AR = Aspect Ratio = (25 /A = 3 


S = Span 
SAT 3 221) Sd m RES 


L D 


bL QA Cus. 
L^ (0/2) VL A D (9/2) VÀ A 


By Figure 62, Ref. (9): 


3 Q 5-4 EM Lb. EI D 
0 0 0 0 0 0. 01 930 
22, 4° 16° 6.40 0. 830 77, 000 0. 096 8, 900 


28.00 — 30" 8.0? 0.980 91,000 0.163 15, 120 





Eos 


el 


ar 


VES > 3 


E OF SHIP 





FIGURE 5 


T, = Side Thrust = L cos (ô - a) - D sin (ô - a) 
D, = Drag = L sin (ô - a)+ D cos (ô - a) 
b - 22.49 
To = 76,000 
D, - 17,400 


5 = 28°: (Breakdown) 


A = 88, 000 
Do = 27, 500 
620 
oO 
= 930 





SYMBOLS: 


E 


Coefficient in series representation of circulation distribution 
Nondimensionalized height of mean line above the nose tail line. 


Slope of the mean line relative to the nose tail line 


Drag coefficient = yr 
TR, 7 oV 


Lift coefficient - lo 
PVS ri /4 A 
Coefficient in series representation of induced radial velocity. 


2 E T 
Thrust coefficient = TO! 


T Rp zP V ^ 
Position vector from vortex element to pivitol point. 
Form Drag 
Induced Drag 
Incomplete elliptic integral of the first kind 
Radial force per unit arc length 
Axial force per unit arc length 
Unit vectors in Cartesian coordinates 
Integral occurring in the series representation of induced radial 
velocity. 
Modulus in elliptic functions 


Complete elliptic integral of the first kind 


2 
u = constant for a particular pivotal point 


Chord length 
Lift force 
Propeller R.P.M. 





2202 


An arbitrary field point where the velocity is to be determined. 
Induced velocity. 

Velocity induced by the bound vortex distribution. 

Radial velocity induced by the propeller. 

Radial velocity induced by the shroud. 

Velocity induced by the trailing vortex distribution. 

Axial velocity induced by the propeller. 

Axial velocity induced by the shroud. 

Radial coordinate. 

Unit vectors in polar coordinates. 

Resultant force acting on the shroud. 

Shaft horsepower. 

Dummy variable = cae Ve 

Thrust deduction. 

Thrust horsepower. 

Side thrust. 

Speed of advance through medium. 

Free stream velocity. 

Functions used in evaluating elliptic integrals. 

vo +t Le = axial velocity of propeller stream. 

Wake fraction. 

Axial coordinate. 

Induced angle of attack at the 1/4 chord point. 

Angle of attack measured from incoming stream to shroud axis in 
plane of symmetry. 

Convergence angle of mean line. 

Circulation per unit length (positive in the positive coordinate 


direction) 








E = 
11/4 


Subscripts: 


B 


n,m,p 


0,0 
1/4 


3/4 


ll 


Circulation (positive in the positive coordinate direction). 


Nondimensionalized axial coordinate, 


Relative rotative efficiency. 

PEOpCliem chicicncy. 

Polar angie. 

Limiting half angle of shroud. 

X/x} ae Nondimensionalized axial coordinate. 
Angular variable of integration = m -~ (9 - 6). 
Axial coordinate. 

Incomplete ellipue integral of menden 
Density of fluid. 

Nondimensionalized radial coordinate 
Amplitude of angular argument in elliptic functions. 


Angular coordinate of the pivotal point. 


Due to bound vortex 


Índices 


Due to propeller 


In radial direction 


Due to shroud 


Due to trailing vortices 


In axial direction 


In tangential direction 


At 1/4 chord point 


At 3/4 chord point 
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APPENDIX A 


DEVELOPMENT OF THE INTEGRAL EQUATION FOR THE 


VELOCITY INDUCED BY THE SHROUD 
\ 







I. Shroud Vortex System: 


e 
(o ^9 x Ue + Sre ^ 
BOUND 
VORTEX 
E d ———— —— M — HQ ee ———— X 





TRAILING VORTEX 


- FIGURE Al. 


Consider ne control volume enclosing the incremental arc of the bound 
vortex distribution lying between the polar angles 8 and 9 + de of Fig. Al. 
Helmholtz vortex theorem states that vortex filaments cannot end in the fluid. 
This requires that the circulation around the vortex filaments leaving the 
control volume must equal that of the vortex filaments entering. This theorem 


implies conservation of vorticity and, in terms of the circulation states that: 


T(@ + A8) + yx(8) 1] /4 A0 = T(0) 





dr, (9) ~ 
10) + (—— ) 20 + 1,40) 1 1 ¡4 A0 - T 9) 


Therefore: 





2:262 


l  dIY(0) 


vx(8) = L2 de 


H. Velocity Induced at a Field Point by the Shroud Vortex System: 


By the law of "Biot-Savart": 


A-2: 


E T (ds XD) 
ir a 


where: dq - incremental velocity induced at a particular field point or 


"pivotal point”. 


ds = directed incremental element of the vortex line. 


9] 759 


- vector displacement from the vortex element to the 


pivotal point. 


Integrating A2, it follows that: 


A*J 
(B) — 
> 1 l'(0)r| ado XD 
qn (1,9, X) = — un] ara 
and 
A-4 


(9 co =) 
arte 2) = ae | | un 
D 


= magnitude of the circulation around the vortex element ds. 


If the pivotal point coincides with a point on the vortex system, one or both of 


the integrals must be evaluated as the Cauchy principle value to exclude the 


singularity at the pivotal point. 
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UI Integral Equations for Velocity Induced bv the Bound Vortex Distribution: 








FIGURE A2 


"P" is an arbitrary pivotal point where the induced velocity is to be deter- 
mined. p has cartesian coordinates (x, y, Z) and polar coordinates (r, $, x). 
By the law of Biot-Savart: 
À.- 4: 
NN de XD 
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In cartesian coordinates: 





EX: 
d9 = da[0i- singj+cosgk] 
T = r [0i + cos $j--sin $ k] 
A-6: 
D = (x)i + [r cos $ - r]/4 COS 01j4 [r sin $ - r]/4 Sin 97k 
do XD= -de | sin 8 (r sin $ - r]/4 Sin 8)- cos O (r cos $ - r]/4 COS 9) ji 
do [x cos 0 )j 
+d9 [x sin 9] k 
EZ 2 LO ul 
|D| = x dur Try 2r 11/4 cos ($ 0) | 
Nondimensionalizing: 
A-7: 
— 29; =A 
T1/4 11/4 
A-8: 
|D| EA (24 o? +1-2pcos (9 E 
A-9: 


de XD = r1/4 (-Locos (2-8) - 1]d9 7 +A cos8d9j+h sino de K y 


Substituting in A-7, A-8, A-9 into A-4: 


I HO) [o cos ($ -9) - 1]d9 = 
* MEIA [+o +1 -2p cos (8-0)? 
dp, = EMI A 


dq. = j 
By ATI y JA [AŽ + p2 + 1 ~- 2p cos Coo 
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E T(0) A sing de 2 
dbz URN 3/2 
4r 11 4 [A^ * p^ *1-2p cos (9 -6)] 
EIN 
Er z das Ed 


> M8) A cos 8 cos $ de 
dg et B 2 2 
4 +1 -2pcos ($ -8)] 


2 


T (8) A sin 9 sin 9 de E 
S oe 
4mxy jg [2-921 -20 cos ($ -9)]/* 


—) 


E = Ar I8) cos ($ - 9) de 
br 4tTr1/4 [A2 -- p? - 1 - 20 cos (8$ -9)] i 


Integrating over the shroud: 





A -10: 
E © 
Ba N: f I'(8) cos (? - 0) de 
br 4T] /4 [A2-- p? 4-1 - 20 COS Co : 
Similarly: 
A-ll: 


2m O 
Zu 1 _ | T(8)L9 cos ($ - 0) - 1]d6 
bx 4nTy/4 J.g [AŻ + p2+1 -2p cos E 2 


de 


diba = dab ` Taj 


» d -L'(8) A cos 0 sin 9 dO 
(0 [nri/4 D * 9^ +1 - 2p cos (8 -0)] 


+1(8) A Sin 9 cos cp de > 
n o e El 
dmx, ld? + (441 - 2p cos (8 -8)] 
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z. _ : AD i I'(8) sin ($ - 0) de 5 
"Y1/4 “e [12+ 02-41-29 cos ($ -9)] 


Note: dn 5 is of little interest in this problem and is included only for the sake 


of completeness. 


IV. Integral Equations for the Velocity Distribution Induced by the Trailing 


Vortex Distribution: 


TA, de 


A 





FIGURE A3 





2012 


By A~]: 
lI  dr«(9 
DD ED m 
E 1/4 
By the law of Biot-Savart: 


>» _ 1 1  dr(8) dg X D 
do = (8 E do 








A-13: 


>. 1 (ar) XD 
D a an DJ? 
D = (x-E)i+(r cos ® ~ r] /4 COS a)j 


+ (r sin $ - r} ¡y sin 6)k 
dé = d£ i 
d£ X D = (0)? - d£ (r sin $ - r4 4 sin 0)j 


+ dí (r cos $ - T] /4 COS 9)K 


172 
|D| = ota - 2r rj/4 cos (3-9) | : 


Nondimensiona lizing: 


Bye /; 











z905 


|D| = ry 74] 6 - e 9 0% + 1 - 2p cos ap? 


-b y 


d—€XD = - r1/4| de (p sin ® - sin 0)j - de (p cos $ - cos 6)K | 


Substituting A-14, A-15, and A-16 into A-13: 





dq+px = 0 

day = darę) A 
m mro d. [Q-e + o? 4-1 - 20 cos (3 -9)] s: 

daz = - __L__ (dT@)\___(p cos 8 - cos 9)de do k 

EA 4TIT ] /4 C de I - cy + o^ + 1 - 2p cos (8-6)? č 


dq. 7 47 izj - dq, ' [cos $j + sin ĉ k] 


en. 1 Are) [sin ($ -6)]dedo r 
m ET ee 
Ir ATA de LA Y t p T 1 - 2p cos (? -0)] s 
dqps = dG, "Ta = dg,.- [-sin 8} +cos $k] 
fee = dr@)  [P - cos (5-8)] de de $ 
Es az zh [( - e) - pg? - 1 - 2p cos (8-9)] $ 


integrating over the shroud: 


- c r y 
Ir ” 4TIE  /4 la o do u a E E In > 
LA - ey t p^-F1-2pcos($-0)] 








230m 


A-18: 
dre” a sis [E ar NO Lee 
-1/4 - [A - €)^ - p^ -- 1-20 cos (8$ - 0)] 
n -19: 
qyx = 0 


Considering A-17: 


ÉL r (0) .. E de 
iTr 7 nri 74 J d 9) imc y 27972 


where k? = $ + 1 - 2p cos (¢-6) can be treated as constant when integrating 
the integral in brackets (ké 20). Integrating the bracketed integral by parts: 
Let u = à - €;du = -de 


LIMITS: e =0 => u=A 


co ER À 
f de E f _ du = u 
O [A - ey + k299/2 A [u2 + 227572 kelu 4 x2172 


NU E \ aa 


: TS 
KAZ +k K 


MEET [. dP(g) sin (8 - 9) de 
-— TÀ S E EO 
Ir ATA -(9 de [o^ +1 -2pcos (5 -6)] m p^ +1 -2pcos(- 91 : 








LI. Lu S 
AT] y4 -0 de [ p? +1 -2pcos ($ -9)] 





aod - 
V. Complete Integral Equations for the Velocity Induced by the Shroud: 
Combining A-10 and A-20: 


dr, o dpr t iTr 





E21: 
» O 
= AE f I(0) cos (3 - 9) de 
Is 2 -0 [A2 + p? 41 - 2p cos (8 - 0)] : 
E oae e o o o 01 n 
mr q Y [(p+1-2pcos(8-0)] [2 2 € 1-20 cos (8 -6)] 
LE P 31) sin(B-e)de 
4TIT 1 JA -0 de [p2 + 1 - 2p cos ($ -0)] 
Ball: 


i: De) Lp cos (6-9) - 1]do 


-— a i 
= q A ———————— 
"xs = dbx 4Tr]/4 Lo [A2 + p2+1 -2pcos ($ - 8) 





Ege 


APPENDIX B 


SOLUTION OF THE INTEGRAL EQUATIONS FOR A PIVOTAL 
POINT COINCIDENT WITH THE TRAILING VORTEX SHEET 


I. Reduction of Eq. (A-21)to Inegrable Form 


For a pivotal point on the trailing vortex sheet, r 2rj/4, Pp” 1, and the 
integrals involving dI'(g )/d9 which arise from the trailing vortices must be 
evaluated as Cauchy Principle Values. 


9 T,., cos (¢-§) de 


S 4TIT] /4 -0 [A2 +2 -2 cos (2 -9)] : 


(3) 
= ar . T 
+ Ar E (9) sin (? - 0) do 


de [1 - cos (3 - 9)] 112+2 -2cos(3 -0)]' 


" ® 
ti 6 Sl) sin($-0)0 
8TIT] /4 E de  [1-cos($-0)] 


Assuming that the circulation distribution can be represented as a com- 


plete Fourier series: 


(9) = ) Ag cos (18) + Bp sin (n9) 
n-o 


The boundary conditions on this representation are: 
a. IY(8) is an even function of 9 for a shroud symmetrical about the 
8 =0 plane. 
b. T(®) = I(-0) = O (8 z riy that is the circulation goes to zero at the 
tips of the shroud. 
By (a) B} = O for all n. 





"368 


r (e) ) An cos (ng) 


n=0 


00 


are) =. ) n Ap sin (18) 


de 
n=O 
R 24: 
2 
Berne A L 


Substituting A-23 and A-24 into A-22 and summing over the integrals: 


n 
-> n=0 > cos (n8) cos (9? - 0) dO 

e A a - 375: 
8427 r1/4 e [4 - cos ($ - o)? é 


n=0 sin (n0) sin (? - 8) d0 
8/2nri y4 -@ L1-cos(8 -0)] [4- cos (8 -o))! : 


00 


r)n 
pou d sin (9-0)d9 _ 
8T r] 4 ^ [1 -cos(é-6)] 


Transformation of variables: 


nð = né - n( - 8) 
cos n0 -» cos n(? - 0)cos nq t sin n(? - 0) sin ndo 


sin nO 7 cos n(? -0) sin nq - sin n(9? - 0) cos no 





Te 


Define: 


since 


v=T-"(%-9); dy = dé 

sin (9-9) = sin y 

cos (5-9) = = COS y 

cos ng = cos n(m- y) cos n? * sin n(rr- v)sin n? 


sin ng 7 cos n(r'- y) sin n? - sin n(Ti- yv) cos n? 


n(T-v) = nT- ny; 

cos (nr) = (-1) 

sin (nr) = (0) 

cos n(m- v) = (-1)! cos nv 4 (0)? sin ny = (-1)! COS ny 
sin n(r - v) = (0)! cos nv - (-1)! sin ny » -(-1)! sin nv 
,.cos n8 -» (-1)![cos nv cos n? - sin ny sin n$] 


sin nO = (-1)! [cos nv sin n? -- sin nv cos n$ ] 


Revisedlimits of integration: 


when 


when 


Note: 


8 =+0; y=T- (8-0) = T+0 -9 = Vo 


0 REUS Vy 


Restriction: vj, v zm 97t8 


Further defining: 


A COS 04, cos (ny) cos y dy 


In = | 


vy,  ([t+cosv] 


8/2 "Ez /4 





A sin (ne) A, 


Cy2 = (-1)" 
8/21 r1/4 
n+l À n sin (n$) A, 
mo A r 
8/2 Tri 
E à n cos (n$) A, 
C4 = (-1) EN 
LES 
4, A Sin (n) A 
Cas = (1) —————— 


STE y4 


(pat n COS (n®)A, 


TTA 


Cnó 


2398 


a V2 sin (ny) cos y dy 
m2 J [2+ cos y] 


vl 


V2 cos (nv) sin v dv 


I, = E cc QUERER 
i V1 [1-cos v][44- cos ] l 2 
me I sin (ny) sin v dv 

n4 — 


w l1-cos v][4 * cos SHE z 


_ ¢ "2 cos (ny) si d 
Ins = | cos (nv) sin v dv 


Y [l+cos y] 


I» sin (ny) sin v dy 
[l+cos v] 


n 


hi6 


= 


Based on the above definitions A-25 reduces to: 


A -26: 
2 8 
dr, 7 r ) ) Com Ium 
n=o m=l 


Solutions to the above integrals are given in the following equations: (n< 3) 





Integral 
lo1 
Hi 
Io] 


31 


192 
112 
152 


132 


lo3 
l13 
l23 
133 


Eq. 
A-43 
A-44 
A-45 


A-46 


N=27 
N 28 
N=29 


A-30 


A-31 
A=32 
A -33 
A-34 
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= 
DR 
an 


Inte gral 


114 
134 


134 


los 
I15 
IT 


135 


loó 
16 
126 
136 


Eq. 
A-47 
A-48 
A-49 


A-50 


A-35 
A-36 
A-37 


A -38 


11:39 
A-40 
A-41 
A-42 





zs 


Note: In evaluating the integrals Ijj and Ij, which involve elliptic functions 


the following relations are used: 


po E ON M 
k DEEP aa 


g= — VAL; 1e k'^ = 
JLF] k 


vn | (5*1) (1 cos vp) 
DES OP > Sin = 2(L+cos vy) — 


Note: See Figure ( A-5) for quadrant selection. 





sn(up;k) = sinpp = sn up; cn up = COS (pp 


dn u = ,/ 1 -k2 sin? ọ 


























E E -Snu _ snu 
ns U = —— tnu = scu E sd u TTE 
E i1 _ enu _cnu 
EU  cnu ,c8u 7 inu EL INI 
ET d _ dnu _ dnu 
ndu = —— ds u- —d deu ==; 
c 2 
V¡(u) = —7 LE(u) - k snucdu] 
k’ 
Va (u) = — [2(2 -k?) E(u) - nr - k^ snu cdu Es Eu +4 -2k?)] 
r 3k’ 


4(2 -x2) Va - 3V1] - k? snu cdu nd*u 


5x1? 


V ¿(u) = 





zt des 


summary of Integral Values: 














E. 2 Im : l 
loi 7 m {u ua "27 [E(uz) E(u2)]] (A-43) 
Iy =I [n 12, x) - mio 2, k) | -—— (uj - u5) - 41911 (A-44) 

J2t+1 1 


l5) = Zr (uy -u9) t 44(1 - 4) LV1(uj) zou NEZ] Ay: LV5(uj) - v5(u5)] 


- 24[h1] - 191] (A-45) 


13} = 449(u - U9) + 1242 - 4) [V¡(uz) - V¡(u))]+ 1241 - LY LV5(u4) - V9(u3)] 


+4(1 - 43% [Vg(uz) - Vg(uo)) 





-24 [1917 -301,1) -2£[191] (A-46) 
mo (A-27) 
COS Vy] 
1 el 
l2 7 2| (t8? 46) ? | (A-28) 
COS V2 
2 3/2 i cos V1 
Dj -4|——rL «4| €*9— -auwe ]} (A-29) 
(L+t)2 a 
= 
cos 
l3; * [2 (82 - 4e 8 - 3) (tH. - 8(9 +22) (ee?) "1 ^ (A-30) 
COS v2 
pc pt = =} COS V] ES 
/ Lt c /L-1 








COS V9 





D 


COS V1 


= 103) 





15 7 2V2+t 
cos Y 








o fK4t 8 N [Y 
cos V2 
cos VI 
133 = 75 | 32 a E 2 En - [193] 
COS Ya 
194 = 0 





2 
l14 = 244+1 E -u +E(u7) = E(u) +355 (snuj cdu, - snu2 cd uz) | 





D, 73 @-1) /TFT {2u - 207 +27 CB(u1) - E(ug)] 
ET [Qa -3 Sy) snuy cdu; - (adu -747 )* 
snug cduz || - 2[14] 
Io, = 8(L-1P 131 lito) - Vq(uj) + | Va(u1) E Va(u2) | 


+254 | va(ug) - Vau) |] - 207-1) U24] - (41) C14] 


l +cos 
E n E 2 | 


lis 7 cos yj ~ cos v - [195] 
Ing = cos? wy - cos? v, - 2115] - [los] 
25 7 COS Ww “COS Ww 15 5 


I35 = ES v cos? "» - 2(125) -3[115] -4L195] 


(A-32) 


(A-33) 


(A-34) 


(A-47) 


(A-48) 


(A-49) 


(A-50) 


(A-35) 


(A-36) 


(A-37) 


(A-38) 





Ba 


Ing = 0 (A -39) 
116 = sin yj -sin w -vi tw (A-40) 
ljg * 2sin vw -2 sin y] +3 (sin 2y] - sin 2w) - (£v) * (tw) (A-41) 


Note: Sign is negative if vp lies in the first or second quadrant and positive 


for third and fourth quadrants. 


I6 7 3 (sin? v, - sin? v) - 20196] - [he] (A -42) 
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I. SOLUTION OF THE GEOMETRIC INTEGRALS, I, IN EQUATION A-26 


Integrals I, and I, are elliptic integrals and are most readily solved in 
terms of Jacobian Elliptic Functions. We will defer this step until last and first 
consider integrals I,9, 1,3» 1,5» and 1,6 which can be reduced to elementary 
forms. 


Using the algebraic transformation: 


cos y=t ; sin v = J1 -t2 


-sin y dy = dt 





dy=- 4 
No 
n cos nv = gnít) sin nv - fq(t) 
0 l 0 
1 t PS) 
2 22-1 2t / (1. - t2) 
3 4t? - 3t 3 /a - t2) 


-4(1 -12%) /(1 -t2) 
= (442 - 1) /(1 -t2) 
Considering 1,2: 
V 


2 . 
IB = | sin nvcos yd y 
V1 [44 cos ns 





COS V1 v(t) t dt 


j (1 2) 2 (441 ]8/2 


cos Vo 


(Note reversal 
of limit with 
sign change) 
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honn-0;f(t-0 
A-27: 192 =0 


Fornzl: Le (EJEA j 


it. t dt 


= + 
COS v9 [t+]? - 
By 119, Ref. (3): 


[x P n/2 de EZ. femen u er 


p? 4n 2En 


with a = 4, b= 1, tn=-3. 


COS Y 
A-28: 52 2| 49/7 qux) > l 





cos v9 


For N=2: £,(t) = 2t/ (1-12) 


M V1 (2 dt 


COS Vo certo’? 


Ing = 2 
Integrating by parts with u = t^, du - 2t dt 


dt -2 
ah, — NE "m 1/28 
(+t) d (4 ot : 


) COS Vj cos V1 
p, -2|-—2 a t dt 
To cos w (t+t)!/? 
COS V2 


By 119, Ref. (3): witha=%t, b=1, tn=-l 
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3/2 COS y] 
-t (t+t) 1/2 
A-29: 1I55 24 „ar? - UL+t) | 
22 
(+t)? 72 3 
COS V2 
Born = 3; f,(t)= (4t? - 1) /1 - {2 
COS y cos V 
Ig = 4 | I Ba _ |o td 
COS V2 DES ‘ COS V5 Catt)” 
Integrating the first by parts with: u= t?, du = 3t? dt 
dv = a V= = 
(t+ t)? "mL 
The first integral equals 
"s 913 COS V] "T VY a 
(L + yr 2 cos V2 7T 2 
COS V9 
By 112, Ref. (3): 
2 2 2.2 
M _ 2(8a AE Ja tbx 
Va+bx 15b 
witha=4,b=l: 
8t — 48 pc 
z t 40 2 2 172 | 
(409172 T E (34 44t + 3t ) (L+t) 
cos Ya 


The second integral is just A-28 .', combining: 
A -30:: 


cos Y 
132 ” " (32 - a1 + 822 - 2) (ar) /? - (8° + 24) (L+ mS i 





cos M3 
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Considering Ir3' 


[? cos nvsinv dv 


In3 7 
i 1 [1+cos vj [t+ cos v] ^ 


QUU mos 
cos v. Lt] [a+ 142 


Forn=0: gu(t) = 1 





i _ T i dt 
= O 
Uo cos va)  [1-*t] (1 + 13) 
By 56, Ref. (3): 
ra. in DVO = kb 
v/u 0 /kb° bb’ fu + Vb" 
with 
MS Sa +bt ea Ib 
u=L+t=atbt; az=tb=1 
k=ab’-a‘b = (2-1) 
A 3l: 





p cante SD 
03  /L-1 Ve PEN 


cos 9 
Forn=1: g,(t)=t 


COS V1 Par 


I1 - 
3 Joos va NEIN [rt] 
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Using v and u as defined in Ip3, t=v- 1 





COS V cos cos 
I. = : (v.c LOHN EE - “1 q 
E COS v9  V/U cos w Ju "cosv, v/u 


Rec = Sitt - Verl COS Y] 
E | esr 


COS V2 


A-32: 


cos Vi 
E - [193 | 


cos v2 





Es 


For n = 2: g, (t) =2t2-1= 2(v - 1) 


gp (t) = 2v2 - 4v+1 


L Op ya - diam v1 a X [77 V oa 
: cos va Yu cos v) Vu “cosvy vvu 


sinceu=4+t=1+t-(1-14=v-(1-2) 
dt = dv 


l+cos va V(£-1)tv 


the first integral = 2 


By 111, Ref. (3): 


RADA A 
2 =z - = Va + bx 
er 3b” 


with: x=v, a=(4-1), b=1 
l rFceos vi 


the first integral = - S (2% -2-v)Vv(L-)+v 
l+cos V9 
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reverting back to t: 











COS y 
>) - 5 (20-3-t) VIFT l 
cos v9 
COS V1 
the second integral = -8/ I +t 
COS V9 
> E || ¡cos V1 
the third integral, from Agg = 1 dy 4-1 
vi-l Vh4t $V 4-1 


OS v9 


Combining: 


: nca cos] 
Dg= (E 4) IF Lo, io a 
Alp ient 641) (cos e 


A - 33: 


ED T cos vV] | 
Ing = Ag 7 4) Eat +|1o3 


cos Ww 





Forn=3: gt) = 4t? - 3t 
- A(v - 1 - 3(v - 1) 
4[v3 - 3v2 + 3v - 1] - 3v +3 


gn(t) 


4v3 - 12v2+ 9v - 1 
By the procedure followed in solving 1) 3' 


: l-+cos MI vr T W SUV 
2 | a V(L-l)+v l+cos Mp) 


COS y cos V1 
+9 te dt 





COS v» Vu "cos m vu 
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By 112, Ref. (3): 


f x? dx A 2(8a2 - 4 abx 4 3b^x 2) 


eS 158 va + bx 


with x = v, a= 4- l, b= l, the first integral yields 


l+cos vi 


ee 2 
g| 85-12 7 6 Dv 313 v^ la v 





l+cos Vo 


reverting back to t: 





3 ) COS Vj 
=) qe |812 - 161+8 - 40 +4 - ste tae + 3 + 32 + 6 | Et 
COS v9 
T2 cos V] 
E Eis +10t - 4 tt + 822 - 202 +15 | Vet 
cos 





The second integral yields, from the development of A-33: 


cos V1 


COS Vo 


84 - 3 -t)/Z£-t 





The third integral yields: 


COS V1 
COS v2 


18V 4 +t 





The fourth integral, from A-30 yields: 





Be EN nz 
/L-1 ftrtt+/t-1 |. " 


combining: 
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1 








ES Ez Be 2 _ 35 
lag = 4137 |3t?+10t - 15t - 4 4t +84" - 208 4304 +15 - 45 42 
ae (0 cos y] 

l-1 /tt + /2-1 cos Vy 
_} 8 2.2 2 15 E 
EE || ota OE 4 Ol SO ELT 
15 4 
/& -1 
en At Yo ] | cos wy 
Il —————————— ae 
/L+t + /R-1 
COS V 
À - 34: 
| COS V 
] 
lag = 5 (3 -5t-440+842+104+=) )/2+t ine 
ES i2 : COS V9 





Considering In5: 


COS y 
| men. | g, (t) dt 
n 


[l+cos v] cos v) eH +t] 


For n = 0: g(t) =] 
fe VI dt COS Ww 


=——= = 4n(1+t) 
cos wy Ll1+t] 





COS v9 


ee 
05 I+cos Vo 
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For n=l: En (€) =t 


a M tdt 
157 COS wo (else 


By 32 Ref. (3) 


| a: E a+bx-a tn (a+bx) | 


(a+bx) 3 b 


win x=t, a=b=1 


cos Vj 
cos Y) 





"- È +t- 4n (1+t) | 
1+cos v] 
lis = COS Vj - COS v5 - in Bow 
I]; = cos v ~ COS v9 - [Ios] 
Forn-2: g (t) 22:12 - 1 


COS y COS Vy] 
me, La | 
E Ll +t] cos m 


dt 
(1+t | 








d É (a+bx)? - 2a (a+bx) + a? 4n (a+bx) | 





withx=t, a=b=l, the first integral of l25 =) 


cos V1 





2 
j| zu" - 20-t) 9 4a +t) | 
COS V9 





2598 


combining with -Igs 


2 cos MI 
ls 72 [2 -t -2454n (1+t)} 





COS W 


l+cos v 4 
» DE) E - J 
Ls = (cos V 7 cos v2) 21 COS V] - COS v2)* ín Ir T | 


A-37: 


2 


Ios = COS V] ^ cos? v9 = DERI = Llos ] 


Forn=3: gnít) = 4t3 - 3t 


cos na 
A. f Me gi s [ t dt 
35 COS y) [1+t] COS Vy [itt] 





By the recurrence relation, Equation 2-1112, Ref. (4): 


f xli dx xh a zul 


(a4 bx) ~ np —b (a+bx) 


withx =t, n= 3, a=b = 1 the first integral of I55 yields: 


M Es 
3 (1+t) 


From loc and lis: 


COS Vi 





las = «e Mioa - 2(1-+t) + 4n(1+1) |) m 


COS V9 
= 3|14t - tn(1 +t) | 
cos v» 
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combining: 


leg e =| (cos y Y - (cos 2» | -2| cos? View cos? v2 | 


l + cos v] 
* (cos v) - (cos vy y + cos i 


A -38: 


135 = 3 | (cos v2) -2| bs | - 3 ns | - 4[Tos | 


Considering lg: 


COS V1 f t) dt 


V 
' = | 2 sin (nv)sinvdv _ | a 
n6 = [1 tet) 


l--cos y 
COS v9 


Forn=0: fn(t)=0 
A-39: 

196 = 0 
Eoi: f(t) = vl-t 


E m vl ACD E m V1 I-t P 


I = 
16 ME (I +t] A 


By 60, Ref. (3): 


pv a k 2 -] [-bu 
J 7a = BY a 7 obs 


wih vs1-t-a'rb't; a'21, b'- -1 








u=l+t=at+bt; a=1, b=1 


k = ab’ - a'b = -1-1 = -2 








I vi BER 


a A 


By the identity 


= [Ji -t2 +2 tan! ib 
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d. Vi 


qe 
cos m 


Ei. i+cos x = tan = = (r= x) 
2* “4 T- cos x 


1 [l+cos w l + cos 
tan” A 
“COS Vp 


A-40: 


l16 = sin V1 


Forn=2: f,(t)=2t 4/1 -t2 


cos V1 
Ing = 2 | 
COS m 
by 161, Ref. (3) 
l xdx 202022 


/a2 - x? 


by 171, Ref. (3): 


[ASS - 2 


with x =t, a= 1 


Ing = 2 i £5 )/1-8 - 2 sin” 


"3 - x) 


- SIN Y) tw 


E AO 
Dien a 


cos v9 /1 - 12 


1 t} cos V] 
COS v2 
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156 = 2(sin v, -~ sin v3) + cos vj sin vı 7 cos wsinw - 
Geo) 
er o a 
A 4l: 


: Ir 
156 = 2 sin w - 2 sin vy ub (sin2 w - sin2. E O EN) 


Note: Sign is negative for Vb in first or second quadrants and positive for third 


and fourth quadrants. 


Forn=3: £,(t)= 442 /1 -12 - /1-t2 


Following the method used in I» 6: 


ES VI j cos v] 
lag = 4 | t2(1-t) de 1 ia 
SA A -d COS v /I+t 





COS VI ¿Sdt raf Vy Zar - 95 V1 /l-T 
COS v -t2 COS Vo A -: COS Ww "PES: 


-4 





By 175, Ref. (3): 
3 a 
AA a5 og.» 97 955 
Er > 22 73 (a^ - x^) a2 fa X 


with a - 1l, x- t, the first integral yields: 


-$ (sin? Si sin? v2) + 4 (sin vj ~ sin w) 
The second integral, from Lg yields 


E 126 | - 4 (sin v - sin v) 
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The third integral = - Ea 


A-42: 


4 . 3 
136 E 3 (sind V2 - sin vi} = 0 EA = EN 
Considering L,, and Ip4: 


These integrals must be solved by use of Elliptic integrals and Jacobian 
Elliptic Functions. 

The solutions for these integrals can be found by transforming them into 
algebraic form as was done for the integrals I,9, In3, Iy5, and1,6- Upon 
transformation it is found that they all are of the form solved in series 235 of 
Ref. (5). In utilizing these solutions, the parameters will have the following 
values: 


4 


a=1,b=-l,c=-, Yp = COS Vp 





Ica =D TE. 2 22 039. 9 a 
Era ep 5 508 ~ ger 
mee - VILIS 1 e = 24 


art /IrT X UR 


Ne (a - cXyp 7 b) 2] (L+1)(1+cos Vp) 
E coop 0m ET sin cos em 3 


up * F(oy.k) - Incomplete elliptic integral of the first kind. 
Ep: K) = E (up) = Incomplete elliptic integral of the second kind. 
T (pp, a2, k) = (up, a?) = Incomplete elliptic integral of the third kind. 


Vy = TT-(8+0); wm) = 1 - (p - 0) 
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Note: Due to the restriction on this group of integrals that 


(azy>b>c) Yp>b = -1 


Eu v9 fT 


Considering I,1: 


i = a COS ny COS y dy 
al") Tor cos 72 
v [4+ cos u 


COS V] gp (6) t dt 
cos w (L+t) A-2 fv 
[^ V1 8, (8) (t für 
nl cos v) (i+t)/1-t /Irt SH 


Emo gQ(t)-] 
B. p^ Vo | : ] B 
01 ~ ur) 7= M Fr 
COS W vet Aet /l1+tt /&tt 
Considering the bracketed term: 


i MEC PAG 
t=(1+t)- 2; [x l DE 


COS Vj dt 
Ini = | 


Ma Vj dt 
“cos v» [t-(-4)] /'1-t /t-(-1) /t- (-X) 
By 235. 00: 
r dt F (o, k) 
Te mun i 
/a-g-bX-c P 


with a = l; bz -1, C= za, the first integral in lol =) 
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2 
en o eo 2 





j 


By 235.01: 


Yp dt 
pe og eer ae "ots E ‚k) 
l (eo Tale) b-c) "P 


with a = 1, b = -1, c = -4 the second integral in lp] > 


E 


——— —- LE(o Kk) - E(po, 
DAT LE(p1»k) - E(p2, k)] 


A-43: 


DES L 
lo. = Fepe (en v2 - rea [Bu EG] 
Forn=1: g,(t)=1 
I a d u | dt 
E COS v9 ero Vi=t /l+e /L+t 


Considering the bracketed term: 








e E [(4+t) - 412 _ (L+t)* - 2UL+t) + 42 
(L+t) (1+t) (L+t 
2 2 
m _ 4 "TE 4 
(rt) = 4+t 24 t TY t USE 


iig] 


comparing with loi 


EE d t dt 
Ii = Eu 5 cene c no c E lo] 
COS Ww /1-t /t- (-1) /t - (-4) 
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By 235. 16: 
2 3 
A DD 
bs t™M dt = om m € p Sn u) E 
b Vla-t)(t-b)(t-c) o (1 -kê sn2u)m 


Wwithim=1, a=1l, b= -1,c= -4 








/a+1 ttl 


The integral in I, ; becomes 


-2 p) (1 - tk? sn? u) ee T2 ie (I - ux? sn? u): E 
Zl "o (l - k2 sn? u) HL ʻo (1 -k* sn? u) 
By 340. 01 with c2 - k^, o2 = 41 


mje 
l-4k snu dn = = [a2 - k2)n(», k?, k) + 4k°u] 
l- k? sn2u k 


- |a - à) (o, k?, K) +4u| 





A-44: 
2 ED) 2 2 RL 
I = “ | 1T1(91,k%,k) - (99, k“, k) OREL 
EE taal ne | Ari ^ Jo 


Where, by 111. 06, Ref. (5) 
| Etk) - (k2sin o cosq/V1 - k? sin? | 
TT (cp, k2, k) = I ixscczumK GG C GG 
k’ 


For n = 2: gn = 2t2 -1 


1 D Qe: dt : 
20 os y CD OE e 





Ge 


Considering the bracketed term: 


2t? -t A Ly RE 
[t+t] rt 


— 26-09 - 661-0) 4 4 e - 243 - (Ltt) +24 
Ut 


2(4+t) - 6UL+t) + Ws elec ei 


= 242 + Mt + 27 - 662 - 655 642 1 - [25-47 


212 - 21: + 2:2 E | 


22-244 e )* ER, 
(+t) (L+t) 
Comparing with I4, and Ip]: 


t^ dt 
EN -0-CDJIU- C201 “P(t | [o 


Again citing 235.16, with m = 2 


jf da - Au EU ue * sn uy du 
¿OT ED E) o (l-kKsn u)? 
of) 0 - 4 sn? uy du 


O (1-k^ sn? we 


By 340.02: witha? = k?, 0 = 2k? 
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2 2 
[Hose Tr = M m = [aktu + 24k (k? - ak yv, + (K? - ay v, | 


By 315: (since a“ = k2, Vj = 1)) 


Vi(u) - => [E(u) - k2 sn u cd u] 
K 


2 


1 2 
VW) = ad [2(2 -k2) E(u) - k*%u] -k? snu cdu UE nd u-c4 E 


4(2 -k2) Vo -3Vi - k? snu cdu nd*u 


9 eee 
3 Sk! 


Note: Vi, V5, V3 are used to evaluate both 15;, 131, and I3,. They are most 


readily applied by evaluating them and substituting in their values. 


2.2 
[ede - n uod 2) Vy(u) + (1-2)? V 34) | 
(1 - a? sn 


Combining 
À - 45: 


li = 24^ (ui = Uy) + 41-2) [Vy (uy) - Vy (ug) | 


+ 21-4" [ Vou) - Vo(ug) | - 22 [111 | - [Tox | 


For n= 3: g_(t) = 4% - 3t 


m I Pe ——— 
3 eos Y) O 
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Considering the bracketed term: 


4° - x^ du - 3) 
TUO E 


2 = [(t+t)- 1]? = Ht) - 244 +0) + 42 


D 
C673. L2 eo - gu SE | 








L+t ttt 
wo > 342 
= 4t° - 444% + (42% - 3) [24+ t) - 21} + de 
3.2. , 62 3 42% - 34% 
_ 2 2 | 
= 4° 24 4242 - - 24 4 2t* 1 [24 vs a 


en { 2. 2. ..209-4 
4t" - 2t 124^ - 24x 21^ - 1- 57 





2 
-4 = 
ICE C 
comparing with Ij5, Ij, and lig 


COS v 3 


= wO E A AA : : 
T s | y Mtro 2, [1 | - 31, ] - 24 19 | 
considering the first integral: 


By 235. 16, with m = 3: 


u u 
4 f l 1 - 22 sn? u)> du -4 | 2 (1 = 22 ene du 

U 3 
o  (1- k^ sn u) O (1 - k^ sn? u) 





Ede 


By 340.03, witha? =k2, a,” = 2K?: 


Eos 
ESA = 2 du d LO ue 30k k* Aa?) v3 a2 - uy v, 
(1 -k^sn^u) 


3 
+ (2 - uy vg 


een - 4k2 sn? u)? 


Sn au = [8 u+ 322 0v, «340-0 v, a- o? vs | 
(1-k^sn^u) 


A -46; 
Igy = 403 (uy up) & 1222 a - 0 [Vy (ay) - Vz (up) 
#124 (1-27 | Vo(uy) - Volug) | 
+4(1-4) [Vaup - Va(up) | -24[ 1, -3[ 1] | -24 Lo, | 
Considering I4 


B sin nv sin v dy 
114 u | 1 
vı [l+cos v][t+cos v2 


| 


NZZ a 
COS v [10t] [44 t 2 


COS y] f(t) dt 


Series 235, Ref. (5) will again be used to solve these integrals. 
For n=0: f,(t) = 0 

A-47: 194 =0 

Forn=1: fp) =%(1 -t 


La n P V] zm - 12 A 


COS vo (+t) Y 
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considering the bracketed term: 


1 -t2 Ve Are 1 


Lnd =, 
= 


1+t 1+t = 





7 
— 
pt 


2 /1-t d 


co RECIEN OUI 


By 235. 07: 


y p" 
P La dt = (a-b) 2 
) - - = g cd“ u du 
| (t-b) (t-c) J 


witha=1,b=-1,c=-4 








Uy B u) 
La = 4 | ca? u du E f cd2 u du 
V&+1 “o Vi+1 0 


By 320. 02: 


i cd u du = = [u - E(u) +k snucdu] 
k 
A-48: 


2 
Iya = De u Eu SES E(uj) * 731 (nui cdu] - snug cduz) | 


For n= 2: £,(t) = 2¢/(1-t2 


COS Vy] 
E- f 2t/ 1 il dt 


dere) et 
cos v ^ Lett E 
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Considering the bracketed term: 


Do etu Tet 


= = 2(1+t) /l-t - 2 vit 


vl+t lt dt 


> 
T 


Comparing with I 14 


Ana AH. 
Bi soos W Vit 7 ES 


COS y] 


By 235. 09: 


Yp = Un 
| o = (ab) (0-0) 8? | sd? u cd2 u du 


With a = Hr b = -l, (C = the integral in l4 =) 


2 2% - 1) GDL W sd? u cd u du - 11 sd? u cd? u du | 


By 361.27: 





fsa? ucdud = > ; - -2K Zu + (1+K’2) E(u) 
3k*k" 
+k? (Kk! 2 nd? u-1-k'2) sn u cd u | 
since: 


2 12 5 
1 — (t1) l+k A 2A K. 
x 4 Á k’ 2 Be y "2 — (2-1) 








TR 





2 
| se? u cd? u du = ere | -2u4+ 74 En 


(ne? u- sU )sn u cd u | 


I 





À.-49; 


Ly = ae 1)v2+1 {2u2 - 2uj + | Eu) - E(w) | 





20 me ZU 
er | (nd ul - ¡TJ snuy cduy 





= (ne? U5 - EE )snug cdu, | | x 2, | 
For n= 3: £,(t) = (42 -1)/1-t2 


"E e Es uri ied dt 


1 
COS v ptt (2 + t)2 


Considering the bracketed term: 


(4t2 - 1) 4/1 -t2- 91 (QtvV1+t2) vi c 
(+t) (1 +t) = 


Substituting 2t = 2(1+t) - 2 =) 


Do R 
(1+t) (1+t) 


Comparing with 14] and 149: 


COS y Vi-t/IRt o. 
34 7 f UR. E /L+t 2m] [ha | 


with t= (4 +1t) - 4, the integral in I¿3 becomes 
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El COS y 
a A A 
u SER) vt-(-4) 


The second integral is the same one as was solved in I,» 


a ie eo ieee ae " DNE 


cos m '  Vt-(-4) 


24] T4 | - 44 4| 


11 


Considering the first integral: 


By 235. 14: 


> u 
P P 
J J- (a -t)(t-b)(t - c) dt = g(a - b) (b - c)? k2 le sd^ucd2und^udu 


with: a= l, b= -l, c = -4, the first integral =) 


úl u2 
4(—* ) 2) A - 1) GI! sd? u cd? u nd^ u du - | sd? u cá? una? udu | 
V L41 a O O 


By 361.18: 


2 


[sd2uca? und udu = AL ptor) I4 - k'2 Ie | 


where Í= V] 
l4 = V2 as defined in I2] 
ee 

A-50: 


I5, = 800 - 12 CET (V3 (3) - V Q3) E | Va(u]) - Va(u)) | + 54) 


| Va(ug) - Vau ] ) - 269 [154 | -(44+1) | Xy q ] 
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APPENDIX C 


COMPARISON OF RESULTS WITH KNOWN RESULT FOR A 
COMPLETE SHROUD 


By the Kutta Joukowski Law: 


(Equation P-2) F(9) = -p Vxp (1,6, 0) Ig) 


Taking the total lift in the 9 = 0 direction: 
[t 
T(0) = -p J Vyp 8, 0) Pig, cos 6 1/4 ds 


In the absence of a propeller 

Vxp = Vz 
TT 

MATO) = DV, 21/4 f Dg) cos 9 de 
“TT 


i23: 
Tg) = 3 An cos (nd) 
n=0 


A 
cos nO cos 0 dé 
TT 


MI) = Son 1/4 2 | 


Because of orthogonality: 


TT 
| cos nð cos 8 dð = On Á 1 
u 


Y 
, TO) = -p Vat, 744, | cos%8 de 
-T 





Sys 


A751: 





T(0) = -p V4 ri/4 Ai" 


To find the side thrust on a complete ring, in the absence of a propeller, 
we need to evaluate Al. 

For a complete shroud, of constant chord, at an angle of attack, the side 
thrust is independent of the section shape (Ref. 6). We will consider the case 


of a complete shroud consisting of a frustrum of a cone as shown in Fig. Al. 





FIGURE As 


The shroud has a chord of L, a radius at the 1/4 chord point of r] /4» a con- 


vergence angle B and is set at an angle of attack, in the 9 =0 plane of a(0) = Ay: 
CO 


ad) = -B+ 0, cos $ 





le 


*s] - rv. [+8 - Ay Cos $) 


EA -26: 
o 6 
ES x > ) Cnm Inm 
n=o m=l 


In evaluating the integrals 1,,, we must abide by the restriction that 


Vy» Y) ÆT 
where 

VI = m- (+0) 

2 T- (9,0) 


Since ® = m, for a complete ring, 


Vi = -9 
D m-t; V9 7 Vj - 2T 
Note: In evaluating the elliptic functions Elo p»k) and F@p,k), the following 


relations are used: 


E(-~,k) = -E(@,k) 
F(-~,k) = -F@,k) 
E(mT +(p,k) = 2m E(k) + E(@,k) 


F(mt +o,k) = 2m K(k) + F@,k) 





je 


hiere: 
E(k) = Complete elliptic integral of the second kind, 
F(k) = Complete elliptic integral of the first kind. 


Determination of ọ] and py from y, and w: 


oat J (L +1) (1+cos vp) 
pee sm 2(£+cos vp) 


2 sin” p (L+cos vp) = (L+1) (l+cos vp) 


= 4+4cos vp + 1 + cos vp 


2 E 
(& cos vp) (2 sin Qp-1)7 1+%cos vp 
l + 4 cos vp 
-cos Qop) = TF cos vp 


when vp changes by 21, | 2cpp| changes by 25. To find the sign, assume 
4=2. 


(1) let VP = 0 


- COS 2 PP. al 


(2) let VP, - m/3 


- COS epp = iti = 0.8 


2% 





Di» 


TT o 
2pp, = 1-37" PP, PP, 18 
(3) let vpg = 0/2 
-cos2pp = 3; cos 2pp = -3 


2pp = 27/3 Pp = n/3 





2738 


(4) let D. p 
4 
-cos 2 — —— 
YP4 1/1 l 


op, = 1 ep, = 0 


when 


vp changes by +2 Tr 


p changes by -n as shown below: 





PIGURE A 9 


E(pp - m, k) - E(pp,k) = -2E(k) + E(pp,k) - E(pp,k) 
E(®p - mk) - E(pp,k) = -2E(k) 
Similarly 


E(pp 7 m K) - F(op, k) 7 -2K(k) 


Since in this case 
va = Zu 
» pn 9-4 -n 
E(u2) - E(uj) = Ep] - mk) - E@1,k) 


E(uz) - E(u]) = -2E(k) 
Similarly 


un > U] E -2K(k) 
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Dy A-43: 
DE 42 K(k) - i E()] 
5 
where 
4242 M 
4 = > 3 Ua, 
2 
sty = 24241 = 
2 (2 +4) 
pu. t-1 
EE e LT 


244 
(+1) = mo N 
Dr 





ZO 
2-1 EC 
lor E cs E ) ElK) + K(k) | 
= 4/2 | - G) Eq) +] 
) Com lom = 901 101 
m=1 


A^ 4-2 AK(k) 
t EL Gg poor vd 
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Bor n= I: 


5 E D 

















By A-44: 
: = 2 2 ce 
24 
: [ +2K(x)] - 4 [101] 
Au! (k)] 101 
Since 
p 4-l- 
R= (Gr) 
li = Ue E(k) + GE y) 50%) - 24 K(k) - 24K(k)) 
2 
St 1) | EG) - 2 K(9)] 
By A-48: 


IJa = 2/2 +1 [+2 K(k) - 2E()] 

4 
I, = —— [-Ü+1)Ek) +(t+1)K(k)] 
14 fT ( ) E(k) + ( ( 


By A-40: 





EE 


) Cimtim * €11l11 * €14h4 * C€16l16 





m= 
h A, cos 0 
=S eee Ty +h, t+ ng ] 
8/2nr, 
1/4 
a ee (HE) 20 - 4-1) K00) 
del 14 JEZI 4-1 
2 
Lra? D Eod. HANC, 
4 = E oU | Le l = — F url ceeias Ir je 32 


Lith. F («eg E(k) - 12 K(x)) 


+4 
A -54: 
6 A,cos $ 3 
E ———— iin T 2) A T 
Cimlim = + E(k) ——— tT 
ES 244 ah? +4 
Similarly, it can be seen from the forms of C,,, and 1,,, that for a com- 
plete shroud: 
6 Ay cos 2 9 
), Com tom 7 7 "rij. NN 
m=1 
and in general: 
DO: 
6 (-1)0+1 A, cos nd 
Be... 1 


Tr 
ET 1/4 
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where H,(\,k) is some function of à and k to be evaluated. Substituting A-53, 


2294, and A-55 into A-26, and equating A-52 to A-26: 


A 
Va [B - dy cos 8] = - 2 | - na (9 *- 
1/4 21/12 +4 
A, cos $ 
„1° (‚7 022? aep. sn -- 
1/4 A 
= is cos nó 
+) ——— [10,6] 
n=2 
Equating coefficients: 
-V_ Tr B 
A = 1/4 
[ M +2 rud 2 AK(K) ae | 
2. J/r2+4 o a 
1 pur uox cc M 
(+ a +2) p 1? K(k) +2} 
244 T 4 
Aa = 0, n22 
Substituting A; into A-51: 


T(0) = maS o A 
(+ TES url 


4 / 2 + 4 


XP 


3 TA 
au 





oru 





EA M Gry 


and utilizing a lift coefficient analogous to that used in Ref. (2) 


E  — —1 
(0/2 V; Lrijg  PVaAr a 
E “N a 
* COELI E 
A-56: 
CL + 4112 


on LE ETA 1% K(k) 
A E(k) - 
{ feta i N2+4 ud 


This result is equivalent to that obtained by Weissinger in Ref. (1) and, as 
shown in Figure A-6 is very close to the result obtained from the more exact 


theory by both Weissinger and Morgan in Ref. (2). 
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